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THE RHEOLOGY AND COMPOSITION OF CRYOVOLCANIC FLOWS ON ICY SATELLITES; 

Jeffrey S. Kargel, U.S. Geological Survey, 2255 Gemini Dr., Flagstaff, AZ 86001 

The rheologic properties of terrestrial lavas have been related to morphologic features of their flows, such 
as levees, banked surfaces, multilobate structures, and compressional folds. These features also have been 
used to determine rheologies and constrain the compositions of extraterrestrial flows. However, with rare 
exceptions, such features are not resolvable in Voyager images of the satellites of outer planets. Often only 
flow length and edge thickness of cryovolcanic flows can be measured reasonably accurately from Voyager 
images. The semiempirical lava-flow model presented below is a renewed effort to extract useful information 
from such measurements. 

The model supposes that a flow moves downslope until its chilled crust attains a critical fraction of flow 
thickness. The solidifying crust is assumed to thicken at a rate proportional to the inverse square root of 
time, analogous to the solidification of the crust on a lava lake. The model does not explicitly account for 
fracturing and foundering of crust, for suspension of crystallized solids in the liquid underlying the crust, nor 
for extrusion of fresh, hot slurry at the leading edges of lava flows. Therefore, this model does not 
accurately portray the behavior of real lava flows. Despite this simplistic approach, the model successfully 
relates the Newtonian viscosities of pure phase liquid lavas to actual flow thicknesses and lengths on Earth 
and the Moon, as shown below. An accurate analytical solution to the problem of the freezing crust (the 
Stefan problem) was presented by [I]. The time, L, for the crust to thicken to the critical fraction of flow 
thickness defines a flow emplacement time. A characteristic flow speed, u, is defined as u = L/t c , where L 
is the length of the flow. No unique critical ratio of crust thickness to flow thickness, T c /T f , has been defined 
by actual measurement or by analytical methods. Setting T c /T f = 1/6 forces the model of one particularly 
well documented lava flow, Kilauea’s Royal Gardens flow of 1983, to reproduce the observed characteristic 
flow speeds and the same effective viscosities as given by [2]. Henceforth, it is assumed that this ratio holds 
for all other flows. The applicability of this assumption to icy flow is supported by laboratory data for 
ammonia-water slurries [3]. Although the exact physical significance of this ratio is unclear, it might relate to 
existing concepts of threshold values such as the Bingham yield strength. 

An effective viscosity (the ‘morphologic viscosity’) is estimated by combining the flow speed calculated in 
the manner just described with the analytical solution presented by [1], relating the viscosity of a fluid flowing 
down an inclined plane to other parameters, including flow speed: 

Vmozp- (Pi^sinaTf p c C) / (288LA, 2 Jc) (1) 

where /* m ^ is morphologic viscosity, p, is the density of liquid, a is surface slope, T f is flow thickness, 
p c is the density of the solidified crust, C is the heat capacity of the crust, L is the length of the lava flow, k 
is die thermal conductivity of the crust, and X is a term given by iteration of the transcendental equation (from 

e- k2 /Xern-Q^/C(T m -T 0 ) (2) 

where Q is the heat of fusion of the lava, T m is its liquidus temperature, and T 0 is the ambient surface 
temperature. All temperature-dependent material properties of the crust are averages over the temperature 
interval T 0 to T m . 

Figure 1 shows the morphologic viscosities of 97 flows plotted against measured or estimated viscosities of 
these lavas at their liquidus temperatures. Morphologic data for these flows and other relevant information 
were obtained from literature sources [2,4-11]. These flows include a wide range of terrestrial silicate flows, 
two sulfur flows, one carbonatite flow, and a lunar flow. They range in length from 5 m to 300 km and in 
thickness from 3.5 cm to 300 m. The least-squares fit has a slope of unity and an intercept of about 10 5 . 

The modeled morphologic viscosity correctly tracks the actual liquidus viscosity (Fig. 1) even though the 
mechanics of real lava flows are far more complex than modeled. The dispersion of data about the 
correlation line indicates a la uncertainty of 1.25 log units in the liquidus viscosity thus determined. A great 
deal of important physics probably accounts for the dispersion, but for present purposes these effects may 
probably be ignored. 

Calculated viscosities of three flows on Triton and one on Ariel are listed in Table 1. Both viscosities 
calculated for the flow on Ariel are several orders of magnitude less than calculated by [12]. The liquidus 
viscosities of Triton’s and Ariel’s flows are in the range exhibited by multi-component ammonia-water liquids 
[3,13]. The liquidus viscosities are also much greater than the viscosities of water (Id 3 Pa-s) and brines (10' 3 
to 10 1 Pa-s), and they are much less than the viscosity of ice (warm ice viscosity ~ 10 3 Pa-s). These 
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viscosity comparisons do not necessarily imply that flows on Ariel and Triton are composed of ammonia- 
water liquids, but such a composition is cosmochemically reasonable. 
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TABLE 1. Calculated flow viscosities. 


Flow name 

Satellite 

Morphologic 

viscosity 

Liquidus 

viscosity 

(Pa-s) 

Possible lava 
compositions 
(Pa-s) 

Unlikely lava 
compositions 

Crater dome flow 
lat 14N long 39 

Triton 

1.7 x 10 10 

1.7 x 10 s 

H20-NH3-GKWE* 

Pure water, brine, 
pure arrincniawater, 
or pure ice 

Lob ate flow 
lat 9S long 57 

Triton 

3.1 x 10* 

3.1 x 10 3 

H20-NH3-GKWE" 

Pure water, brine, 
pure arrmxm-watr, 
or pure ice 

Lake spill flow 
lat 20N long 20 

Triton 

1.8 x 10 7 

1.8 x 10 J 

H20-NH3 or 
H20-NH3-GKWE* 

Pure water, brine, 
or pure ice 

Kewpie-Brownie 

flow** 

Ariel 

2.3 x 10 9 
or 4.6 x 10 s 

2.4 x 10 4 
49 

H 2 0-NH 3 -GKWE* 

H 2 0-NH 3 or 
HoO-NFL-GKWE* 

Pure water, brine, 
or pure ice 


* GKWE = God-Knows-What-Else. Might include methanol, ammonium sulfide, sodium chloride, or other substances 
providing substantial added freezing-point depression and increased viscosity [3, 4], 

♦♦First entry based on assumption of axial eruption; second entry, on assumption of longitudinal flow. 



Figure 1. Calculated morphologic viscosities and estimated liquidus viscosities of 97 lava flows. 



